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Hurricanes cause severe impacts on forest ecosystems in the United States. These events can substantially 
alter the carbon biogeochemical cycle at local to regional scales. We selected all tropical storms and more 
severe events that made U.S. landfall between 1900 and 2011 and used hurricane best track database, a 
meteorological model (HURRECON), National Land Cover Database (NLCD), U. S. Department of 
Agirculture Forest Service biomass dataset, and pre- and post-MODIS data to quantify individual event and 
annual biomass mortality. Our estimates show an average of 18.2 TgC/yr of live biomass mortality for 1900- 
2011 in the US with strong spatial and inter-annual variability. Results show Hurricane Camille in 1969 
caused the highest aboveground biomass mortality with 59.5 TgC. Similarly 1954 had the highest annual 
mortality with 68.4 TgC attributed to landfalling hurricanes. The results presented are deemed useful to 
further investigate historical events, and the methods outlined are potentially beneficial to quantify biomass 
loss in future events. 

Hurricanes in the North Atlantic, are major natural disturbances that are part of life in the southern and 
eastern United States causing severe impacts on forest ecosystems in the region 13 . The intensity of impact 
to forest stands is often highly variable, ranging from short-term to long-term. Short-term impacts may 
include defoliation and branch break or total blowdowns. Long-term impacts may include rapid or slowed tree 
mortality, changes in successional direction, increased species turnover and age class diversity, and faster biomass 
allocation and carbon storage 1,4-6 . Alterations to forest composition and structure, and transfer from live trees to 
dead carbon pools are certain regardless of the inter- and intra-annual frequency and intensity of hurricanes 2 ' 7-9 . 
Several studies have suggested that the frequency of hurricanes would decrease but their intensity may increase in 
some locations as global climate warms in the future 810 . The increased hurricane intensity would certainly destroy 
more forest biomass whenever they pass over the forestland. 

Intense hurricanes can cause substantial damage to a forest but their size alone is not enough to determine how 
much live tree biomass results in as debris. Forest types, wind speed, and heavy rainfall followed by flooding are 
the key factors for most severe storm damage 11 to forested areas. Evergreen forests are generally more susceptible 
to hurricane damage than deciduous forests. Susceptibility of forest type to hurricane-related damages also varies 
for different types of damage, including breakage, and uprooting 1214 . In addition, trees planted outside their 
natural range may be more susceptible to damage 15 . Immediate damage to hardwood forests range from defo- 
liation to total blow down; however, many species may recover quickly with releafing and sprouting often 
occurring within 3-4 weeks 16 . 

The role of individual tree mortality through its disturbance on regional carbon balance is complicated by the 
strong spatial and temporal nature of slow-in and rapid-out of carbon flux 17,18 . Thus, it is important to understand 
the spatiotemporal variation of disturbance and its effects on forests and carbon dynamics to quantify current 
mortality and to predict future change 17 . Given the variations in direction, duration, and distance of wind speed 
relative to the eye of the storm, predicting forest mortality from hurricane events is difficult. Different responsive 
natures of different tree species, age classes, soil profiles, as well as topography add to the difficulty 3 . 

As the North Atlantic basin has experienced an increase in hurricane activities with some intense events in the 
recent past, various efforts were undertaken to quantify carbon dynamics from these activities. Some studies were 
primarily targeted to study biomass damage from individual events, whereas others were involved on extended 
spatial and temporal scales 2-6,19-20 . Chambers, et al. 19 used Landsat data, field-based information, and an empirical 
model to quantify biomass loss from Hurricane Katrina. Zeng, et al. 2 applied similar approaches with the addition 
of simulated wind field for each hurricane event to model biomass loss and carbon release from 1851 to 2000 for 
all US landfalling hurricanes. Nielsen 6 used Moderate Resolution Imaging Spectroradiometer (MODIS) data, and 
environmental and meteorological information to predict the level of damage to forests in the Katrina impact 
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Figure 1 | Forest areas were derived from NLCD 2001. Hurricane best tracks and category data are from HURDAT database for 1900-2011. Categories 
are based on the Saffir-Simpson scale and TS stands for tropical storm. Best tracks were cut off after they decreased below a sustained weed speed 
of 18 mph over land. Map generated using ArcGIS 10.1. 



Table 1 | Ten most destructive hurricane events with dates and categories at landfalling. Numbers in parentheses show minimum and 



maximum range 


Name 


Landfalling Date 


Category at Landfalling* 


Biomass mortality (TgC) 


Camille 


8/15/1969 


5 


59.49 (41.42 ~ 77.58) 


Donna 


9/10/1960 


4 


51.48(35.80-67.16) 


Hazel 


10/15/1954 


4 


47.39(32.27-62.51) 


Okeechobee 


9/17/1928 


4 


41.22 (28.78 - 53.69) 


Elena 


9/1/1985 


3 


38.42 (26.98 -49.86) 


Katrina 


8/29/2005 


3 


36.03 (27.75 -46.29) 


Gracie 


9/29/1959 


3 


33.85 (23.47-44.21) 


Diana 


9/13/1984 


3 


33.43 (23.18 -43.68) 


Hugo 


9/22/1989 


4 


30.47(20.91 -40.05) 


Frederic 


9/13/1979 


3 


30.4 (21.22 - 39.58) 


*based on Saffir-Simpson scale. 
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Figure 2 | Histogram showing biomass mortality distribution by years (A) and by events (B) between 1900 and 2011. 



Table 2 | Ten most destructive years and number of hurricane 
events with biomass mortality values. Numbers in parentheses show 
the minimum and maximum range 

Year Number of Events Category >= 3 Biomass mortality (TgC) 



1954 


4 


3 


68.37(46.57-90.18) 


1960 


4 


1 


65.95 (45.98 ~ 85.93) 


1969 


3 


1 


63.76 (44.38 ~ 83.18) 


1985 


6 


1 


62.66 (43.94 ~ 81.38) 


1916 


9 


1 


60.59 (42.18 ~ 78.99) 


1998 


5 


0 


59.17(41.23 -77.11) 


2005 


8 


4 


56.42 (41.56-72.61) 


1979 


5 


1 


55.86 (38.92 - 72.8) 


2004 


8 


3 


46.24 (31.52 -58.81) 


1906 


4 


1 


44.59 (31.23 -57.96) 



region in the states of Louisiana and Mississippi. McNulty 11 com- 
bined historical and current forest damage and hurricane informa- 
tion, and analyzed loss of forest biomass through hurricane events 
using timber loss vis-a-vis the proportion to leaves, roots, and stems. 
After Hurricane Andrew, various studies were conducted to evaluate 
biomass damage in different forest stands from Florida to the Gulf 
coast; however, none of these studies focused specifically on quan- 
tifying carbon from the biomass loss 16,21,22 . These studies recognize 
the important role that hurricanes play in forest biomass mortality 
and carbon cycles. However, none of this previous work proposed a 
simple approach to quantify biomass loss from individual hurricane 
events or produced a live biomass mortality dataset for large number 
of historic hurricane events. We propose a simple method to quantify 
live biomass mortality datasets for individual hurricane events as well 
as annual mortality from hurricanes using hurricane best-track data, 
and the National Land Cover Database (NLCD) as inputs. 




1900 1920 1940 T9f5rj~ 1980 2000 

Year 

Figure 3 | Range of annual biomass mortality (light grey shade) and range of the 10 -year average (dark grey shade) between 1900 and 2011. The mean 
of the 10 -year average is the black solid line, the mean of the 10 -year average adjusted with forest area is the red solid line, the mean of the 10 -year average 
adjusted with biomass stock is the green line, and the mean of the 10 -year average adjusted with a combination of forest area and biomass stock 
and FIA estimates is the red dotted line. 
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Figure 4 | Long-term average biomass mortality from hurricane events for CONUS over different time periods. Map generated using ArcGIS 10.1. 
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Table 3 Annual average biomass mortality divided into 25-year 
periods plus the last 1 2 years of the study period. Numbers in 
parentheses show the range of minimum and maximum values 

Time Biomass Mortality (TgC/yr) 


1900-1924 


18.49 (12.93 ~ 24.05) 


1925-1949 


15.04 (10.63 ~ 19.30) 


1950-1974 


22.61 (15.72 ~ 29.49) 


1975-1999 


17.75(12.37-23.12) 


2000-201 1 


15.46 (10.59 ~ 19.90) 


1900-201 1 


18.15(12.66~23.5j 



This study covers hurricane events that originated in the tropical 
region of the Atlantic Ocean and made landfall in the conterminous 
United States (CONUS). The US hurricane impact zone includes the 
CONUS coastal region, ranging from the Gulf of Mexico to the 
Florida peninsula and along the Atlantic coast. Twenty-five states 
in the region, including the Washington DC area, have experienced 
impacts from historical hurricane events (Figure 1). 

Results 

Biomass mortality by individual hurricane events. Hurricane 
Camille, classified as a category 5 hurricane on the Saffir-Simpson 
scale that struck the coasts of Louisiana and Mississippi in 1969, 
causing 59.49 TgC in mortality of aboveground live biomass 
(AGB), the highest amount in US history. The other most 
destructive hurricane were Donna in 1960 with 51.48 TgC, Hazel 
in 1954 with 47.39 TgC, and Okeechobee in 1928 with 41.22 TgC 
(informally named after the event 23 ). These events were classified as 
category 4 hurricanes when they made landfall (Table 1). These 
hurricanes caused high biomass mortality because they were 
passing over largely forested areas. Only 10% of the 330 hurricane 
events in the study period caused £30 TgC biomass mortality, 
whereas, nearly 60% of these hurricane events destroyed £5 TgC 
(Figure 2A). 

Annual biomass mortality from hurricanes. Between 1900 and 
2011, the average annual loss of AGB was calculated as 18.2 TgC 
(the range was between 12.7 and 23.5 TgC). The highest number of 
hurricanes making landfall in the United States that were greater 
than 35 mph occurred in 1916 (nine), followed by 2004 and 2005 
(eight each year); however, the highest annual AGB loss was highest 
in 1954 with 68.37 TgC. The three years with the highest number of 
landfall hurricanes (1916, 2004, and 2005) included three or more 
hurricanes rated as category 4 or higher. In 1954, only four 
hurricanes made landfall in the U.S.; three of them were rated as 
category 3 or higher. The other most destructive years in terms of 
AGB mortality were 1960 and 1969 with 65.95 TgC and 63.76 TgC, 
respectively (Table 2). In 1922, no hurricanes made landfall on U.S. 
mainland, and there were 12 years (1905, 1925, 1931, 1951, 1962, 
1973, 1978, 1980, 1987, 1990, 2006, and 2010) when <1 TgC of AGB 
mortality occurred from landfalling hurricanes (Figure 3). 
Landfalling hurricanes accounted for more than 55 TgC/yr only 
during 8 years (Table 2). On the other hand, 89 years experienced 
less than 30 TgC/yr AGB mortality annually, while 33 years had 
5 TgC/yr or less biomass destruction (Figure 2B). 

Spatial pattern of biomass mortality. The Gulf coast (coastal region 
of Louisiana, Mississippi, Alabama, and northwestern Florida) and 
the mid- Atlantic coast (Virginia and North Carolina coastal areas) 
(Figure 1), lost higher biomass than any other areas on average per 
square meter over the studied time period (Figure 4), which coincides 
with weather- related forest mortality maps from Forest Inventory 
and Analysis (FIA) inventories 24 . Spatially, the 1925-1950 time 
period had less biomass mortality rate per square meter, even 
though hurricanes making landfalling during that period were not 



far from the long-term annual average (2.94 events per year) 25-26 . On 
the other hand, the periods from 1950-1975 and 1975-2000 had 
below annual average landfalling hurricane events (2.78 and 2.52, 
respectively, versus 2.94 events per year) 25,26 . However, these periods 
had above average biomass mortality rate (see Figure 3 and Figure 4) 
because of some highly intense hurricanes (such as Hazel 1954, 
Donna 1960, Camille 1969, Hugo 1989, and Elena 1985) making 
landfall over largely forested regions 3 . 

Discussion 

The long-term annual average biomass mortality from this study 
(18.2 TgC/yr for 1900-2011) was comparable with the results from 
Zeng, et al. 2 , who estimated an average of 19.5 TgC/yr live biomass 
loss for the 1900-2000 period after synthesizing field measure- 
ments, performing satellite image analysis, and using some com- 
plex empirical models. Annual biomass mortality in this study was 
an aggregation of mortality from all events of the particular year 
on a pixel-by-pixel basis, and biomass mortality from individual 
events was estimated using some simple datasets (forest area and 
biomass stock) and field-measured coefficients. For example, if any 
location was experienced hurricane events more than once in a 
particular year, the maximum destruction from those events was 
taken into account rather than the sum of destruction from all 
events. As McNulty 11 explains, this approach is more reasonable 
because a substantial portion of biomass destroyed by a first hur- 
ricane event would recover by the second hurricane event to hit the 
same area, even if the event occurred within a couple of months. 

Annual average biomass loss was different for each portion of the 
studied time period. The largest loss was estimated for 1950-1974 
with 22.6 TgC/yr (Table 3); this period included the top three most 
destructive hurricane events (Table 1). Overall, 2.78 hurricane events 
per year made landfall during the period from 1950-1974, compared 
to 2.94 events per year for the entire study period (1900-201 1), but a 
higher number of these events were more intense (category 3 and 
above). Most intense hurricanes also passed through highly forested 
regions damaging a large amount of biomass because biomass mor- 
tality is dependent not only on the intensity of hurricanes but also the 
presence of AGB 1 1 . Annual average AGB mortality for the 2000-2011 
period was substantially lower than the long-term annual average 
(15.46 versus 18.15 TgC/yr) but our estimates were similar to FIA 
weather-related biomass mortality data (15.5 TgC/yr) for the hur- 
ricane impact area (Table 3 and Figure 3). The FIA estimation should 
have a slightly higher mortality than this study because it includes all 
types of weather-related forest mortality (ice, wind, hailstorm, tor- 
nadoes, windstorms, and hurricanes) in a single discrete class. FIA 
measurements may have underestimated the actual loss because for- 
est regeneration and succession are rapid after hurricane-like distur- 
bances 11,17 . Additionally, there is temporal scale variability of FIA 
data collection. FIA compiles mortality data through plot inventory 
taken every 5 to 7 years for most states, and to qualify for any given 
disturbance, at least 25% of the trees should sustain damage or mor- 
tality in an acre of land since the previous visit 17 . Particularly, some 
states (for example, Mississippi) have participated in special acceler- 
ated inventories. Following Hurricane Katrina, the entire state of 
Mississippi was inventoried within a 687-day time span 27 since the 
timing of the damage measurement after hurricane events is 
sensitive 3 ' 28 . 

For a long-term study such as this, the dataset availability for 
frequent time intervals is important because a lack of historical spa- 
tial datasets can add a significant level of uncertainty in the results. 
We analyzed the uncertainty using historical forest area and biomass 
storage data reported by Smith, et al. 29 and Birdsey and Heath 30 . Our 
results appear to overestimate by up to 10%for certain years because 
of a lack of historical spatial forest area dataset (red line in Figure 3). 
Similarly, lack of historical spatial biomass stock datasets may have 
offset the result by up to 50% (green line in Figure 3). The combina- 
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Figure 5 | Examples of simulated severity levels and biomass mortality for three events; a, c, and e are damage severity zones of Katrina, Rita, and Ivan, 
respectively; b, d, and f are biomass mortality from Katrina, Rita, and Ivan, respectively. Map generated using ArcGIS 10.1. 
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Stanturf, et al. 3 


Katrina 2005 
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22.5 


40 


Meeker, et al. 38 


Katrina 2005 


1.17 


1.36 


10.47 


30.91 


Oswalt and Oswalt 51 


Rita 2005 


3 


15 


50 


60 


Texas Forest Service 46 


Ike 2008 


3 


10 


25 




Texas Forest Service 47 


Average 


4.03 


1 1.81 


31.49 


49.58 




STD of Means 


1.45 


2.44 


6.01 


7.58 





tion of not having both of these datasets reveals an even greater effect 
(see the blue line in Figure 3). Therefore, at the very least, the avail- 
ability of decadal spatial datasets of forest area and biomass storage 
availability may render lower uncertainty in our results. 

Because the immediate effect of hurricanes in forests is conversion 
of live biomass to dead carbon pool, finding damaged areas after each 
hurricane event is critical to quantifying live biomass mortality to 
illustrate the distribution of biomass damage in the impact zone 5,31 . 
Mapped damage severity zones of this study for some recent events 
were compared with damage severity zones defined by the USDA 
Forest Service 32 (Figure 2). The Forest Service established damage 
severity zones based on field observations of various forest sample 
plots. Almost always, a hurricane caused more damage on the right 
side of its track than the left side (see left panel of Figure 2). As Chen, 
et al. 33 stated, wind and stress could be as much as 25% higher on the 
right side of a hurricane track than on the left side. However, the 
damaging distance from the track could vary from one event to 
another as every hurricane has a different level of wind strength, 
forward speed, and geographic size 34 , as well as various components 
that could affect the hurricane's course and intensity 33,35 . This study 
also found that every hurricane event had a different footprint of 
damage area and level of biomass mortality but the right side of the 
tracks consistently had more impact (see some examples in Figure 5). 

Hurricane Katrina was one of the most widely studied hurricane 
events that landed in the United States. FIA mapped a higher severe 
zone on the right side of the track 32 of Katrina, which coincides with 
this study (Figure 5a). This study also found most of the forest 
biomass damage occurred on the right side of the track 
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(Figure 5b). Our study estimated that total biomass destruction by 
Hurricane Katrina was 36 TgC (ranging from 27.8 to 46.3 TgC), 
which is comparable to the results of Negron- Juarez, et al. 7 , who 
estimated 43.9 TgC (±8.4 TgC) using a relationship between field- 
measured tree mortality and Landsat data. On the other hand, the 
spatial map developed by Chambers, et al. 19 was comparable to this 
study (with Figure 5b) but their estimation of total biomass loss from 
Hurricane Katrina was 102 TgC. Chambers, et al. 19 included surface 
litter, whereas this study and the study by Negron-Juarez, et al. 7 did 
not, and yet, their results appear to overestimate the total biomass 
loss because it is unlikely that surface litter accounts for more than 
twice the aboveground live forest biomass 36,37 . 

The estimated impacts of hurricanes on carbon dynamics are 
highly dependent on the linkage between remotely sensed signals 
and ground observation of tree damage (see Table 3). Often different 
field-based studies suggest different damage levels even for the same 
event. For example, Nielsen 6 reported that Hurricane Katrina's 
severity of forest damage ranged from 15 to more than 50%. 
Stanturf, et al. 3 pointed out that Hurricane Katrina destroyed 
approximately 67% of the trees in the severely damaged zones and 
33-66% and 3-33% in moderately, and lightly damaged zones, 
respectively. An inventory performed by the Forest Service in the 
DeSoto National Forest after Hurricane Katrina shows up to 83% tree 
damage per acre in severely damaged plots 38 . The forest damage and 
severity zones are also influenced by tree species. For example, over 
80% of all pines were overtly damaged by Hurricane Andrew in the 
Everglades of southern Florida 16,21 . After Hurricane Andrew, about 
85% of hardwood hammocks were in some way affected but average 
tree mortality was only 11.5% when assessed after 4 months follow- 
ing the event in the Long Pine Key, Everglades National Park 16,22 . 
Conducing additional field studies to systematically cover more 
events and locations should improve the estimates on the biomass 
loss percentage for each severity class and reduce overall uncertainty. 

Although our approach provides spatially explicit estimates of 
biomass conversion from live to dead carbon pool after storms, it 
has some drawbacks. Uncertainty in the biomass loss such as sal- 
vaging after hurricane events was not considered in this study but 
that could make a significant difference in the values presented. As 
Stanturf, et al. 3 noted, 37% of damaged biomass was salvaged in 
South Carolina following Hurricane Hugo in 1989. Our approach 
can only quantify immediate biomass mortality, yet many tree spe- 
cies tend to die many months and even years after hurricane 
events 311 . In addition, we did not consider the impacts of climate 
and land cover change on carbon pools in the reconstruction of the 
long-term historical carbon loss, which might have overestimated the 
carbon loss especially in the early time periods when forests were 
young and recovering from agricultural abandonment 39,40 . 

Methods 

Characterization of hurricanes. The best tracks for individual hurricanes that had 
wind speeds of at least 35 mph at landfall were retrieved from the National Hurricane 
Center (NHC) data archive 41 . NHC maintains hurricane track datasets beginning in 
1851 with latitude, longitude, minimum pressure, and maximum wind speed at the 
center of circulation documented every 6 hours. The best track information is used to 
simulate wind field with 5-km pixel size using the HURRECON model. HURRECON 
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Figure 7 | Flowchart showing methodology applied in this study; the left panel depicts the remote sensing approach and the right panel shows the forest 
type damage fraction approach. Forest types 41, 42, 43, and 90 are deciduous, evergreen, mixed, and woody wetland respectively; forest type for woody 
wetlands for NLCD 1992 is 91. 



is a meteorological model that can simulate wind field for a user-defmed resolution 
using maximum wind speed, hurricane radius, and geographical locations of the 
hurricane's eyes. The hurricane track data must contain maximum wind speeds and 
geographical locations of the hurricane's centers 1,2 ' 4 whereas hurricane radius is 
calculated using the wind speed and geographical location of the hurricane's 
eyes 42 . 

Remote sensing of storm damage. The MODIS Reflectance product MCD43A4, 
which is adjusted to nadir using a bidirectional reflectance distribution function 
(BRDF), contains visible and infrared surface reflectance at 500 m resolution for each 
16-day period. This product is generated from both Terra and Aqua data to provide 
the highest probability for quality input data. Atmospheric effects have been removed 
from this product as it would have been measured at ground level, so this is more 
consistent for observing vegetation status than other MODIS reflectance 
products 31,43 ' 44 . The MODIS Reflectance products were used to derive Normalized 
Difference of Infrared Index (NDII) because of its sensitivity to vegetation water 



content or weight of water per unit area 31,45 and better estimation of vegetation 
disturbance after hurricanes 31 . NDII is calculated 



NDII = 



(NIR-SWIR) 
(NIR + SWIR) 1 



(1) 



where shortwave infrared {SWIR) at 1.6-2.13 urn wavelength and near infrared (NIR) 
at 0.86 um wavelength were used. MODIS products have two bands with SWIR 
spectra, however we adopted the 2.13 jim band due to the large amount of missing 
observations caused by a serious striping issue in the MODIS Aqua 1.65 um band 31 . 

To estimate biomass damage, we calculated differences of pre- and post-NDII of 
hurricane events as: 



ANDII -preNDII - postNDII 



(2) 



We calculated average and standard deviation of ANDII for the land area that sur- 
passed a wind speed of 35 mph - values that we could use in subsequent calculations. 
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Calculating the impacts of storms on tree mortality. This study followed the FIA 
defined naming convention for forest damage zones to compute severity classes. FIA 
conducts ground inventory to estimate forest damage for each plot as the percentage 
of basal area that suffered no damage, scattered light (Level 1), light (Level 2), 
moderate (Level 3) or severe damage (Level 4 ) 31,46,47 . The percentage of biomass loss 
for each severity class in this study was developed by summarizing results from 
various studies (Table 4). We applied bootstrap sampling 48 with 10 iterations to 
calculate the mean and standard deviation of means in order to quantify the range of 
biomass loss for each event. 

We employed a remote sensing approach (i.e., using MODIS) to estimate the range 
of biomass damage from 2000 to 20 1 1 using average and standard deviation values of 
damage fraction presented in Table 4. Because MODIS data were not available before 
2000, we had to develop a reliable damage fraction for different forest types. As 
Stanturf, et al. 3 and Vanderwel, et al. 17 noted, different forest species and types 
respond differently to various kinds of disturbances. We used individual damage 
severity maps developed through a remote sensing approach from 2000-2011, and 
extracted the biomass damage for each forest type for individual events. We further 
calculated percentage of biomass damage by forest type to establish a relationship 
between forest types and severity level of biomass damage (Figure 6). The percentage 
of biomass damage by forest type from individual events was then averaged and 
applied for quantifying biomass damage for years between 1900 and 1999 {see 
Figure 7); we named this approach forest type damage fraction. 

Land cover and biomass storage. NLCD inputs for 1992, 2001 and 2006 were 
retrieved from http://www.mrlc.gov/fmddata.php. We masked out forestland from 
the NLCD datasets for the hurricane impact region in the CONUS. NLCD 1992 was 
used to derive forest pixels prior to 1995, NLCD 2001 was used for the period 1996- 
2005, and NLCD 2006 was used for the post-2006 time period. Woody wetland was 
also included in this study as forested land with three forest types: deciduous, 
evergreen, and mixed. 

CONUS forest live biomass data were downloaded from http://fsgeodata.fs.fed.us/ 
rastergateway/biomass. The forest live biomass dataset was validated with FIA field- 
measured data for randomly selected plots from 65 CONUS mapping zones. The 
correlation coefficient of the observed field- measured data and the live biomass 
datasets ranged between 0.79 and 0.31 among mapping zones 49 . The unit of the 
dataset is converted to carbon equivalent by multiplying by 0.5 2 ' 50 . 

Uncertainty analysis. Given the lack of historical spatial datasets to map forest cover 
and biomass prior to 1992, we postulate that a propagation of uncertainty was 
introduced. To map the uncertainty, we took historical total forest area and biomass 
stock, and calculated the ratio of changes from the base year of 2000 reported by 
Smith, et al. 29 and Birdsey and Heath 30. The periodic coefficient values were then 
interpolated to develop annual values, which were then multiplied to annual biomass 
mortality to derive adjusted values. 

We assume that coefficient values for severity classes also introduced some 
uncertainties as these values were derived from smaller pool of literature and no 
climatic variables were considered in this study. The results would have been different 
if more literature values were available and climatic variables were incorporated, 
however, this is work for future research. 

1. Boose, E. R., Foster, D. R. & Fluet, M. Hurricane Impacts to tropical and 
Temperate forest landscapes. Ecol. Monogr. 64, 369-400 (1994). 

2. Zeng, H. et al. Impacts of tropical cyclones on U.S. forest tree mortality and carbon 
flux from 1851 to 2000. Proc Natl Acad Sci USA 106, 7888-7892, doi:10.1073/ 
pnas.0808914106 (2009). 

3. Stanturf, J. A., Goodrick, S. L. & Outcalt, K. W. Disturbance and coastal forests: A 
strategic approach to forest management in hurricane impact zones. For. Ecol. 
Manage. 250, 119-135, doi:10.1016/j.foreco.2007.03.015 (2007). 

4. Boose, E. R., Serrano, M. I. & Foster, D. R. Landscape and regional impacts of 
hurricanes in puerto rico. Ecol Monogr. 74, 335-352 (2004). 

5. Lugo, A. E. Effects and outcomes of Caribbean hurricanes in a climate change 
scenario. Sci. Total Environ. 262, 243-251, doi:Doi 10.1016/S0048- 
9697(00)00526-X (2000). 

6. Nielsen, E. M. in Eighth Annual Forest Inventory and Analysis Symposium. 
307-316 (2009). 

7. Negron-Juarez, R., Baker, D. B., Zeng, H., Henkel, T. K. & Chambers, J. Q. 
Assessing hurricane -induced tree mortality in U.S. Gulf Coast forest ecosystems. 
/. Geophys. Res. 115, doi:10.1029/2009jg001221 (2010). 

8. Bender, M. A. et al. Modeled impact of anthropogenic warming on the frequency 
of intense Atlantic hurricanes. Science 327, 454-458, doi: 10. 1 126/science.l 180568 
(2010). 

9. Goldenberg, S. B., Landsea, C. W., Mestas-Nunez, A. M. & Gray, W. M. The recent 
increase in Atlantic hurricane activity: causes and implications. Science 293, 
474-479, doi:10.1126/science.l060040 (2001). 

10. Emanuel, K., Sundararajan, R. & Williams, J. Hurricanes and Global Warming: 
Results from Downscaling IPCC AR4 Simulations. Bull. Am. Meteorol. Soc. 89, 
347-367, doi:10.1175/bams-89-3-347 (2008). 

11. McNulty, S. Hurricane impacts on US forest carbon sequestration. Environ. 
Pollut. 16, 817-824 (2002). 

12. Foster, D. R. & Boose, E. R. Patterns of Forest Damage Resulting from 
Catastrophic Wind in Central New England, USA. /. Ecol. 80, 79-98 (1992). 



13. Jalkanen, A. & Mattila, U. Logistic regression models for wind and snow damage 
in northern Finland based on the National Forest Inventory data. For. Ecol. 
Manage. 135, 315-330 (2000). 

14. Peterson, C. J. Catastrophic wind damage to North American forests and the 
potential impact of climate change. 5c/. Total Environ. 262, 287-311 (2000). 

15. Foster, D. R. Species and Stand Response to Catastrophic Wind in Central New 
England, U.S.A. /. Ecol. 76, 135-151 (1988). 

16. Armentano, T. V., Doren, R. F., Piatt, W. J. & Mullins, T. Effects of Hurricane 
Andrew on Coastal and Interior Forests of Southern Florida: Overviewand 
Synthesis. /. Coast. Res. Special Issue No. 21, 111-144 (1995). 

17. Vanderwel, M. C, Coomes, D. A. & Purves, D. W. Quantifying variation in forest 
disturbance, and its effects on aboveground biomass dynamics, across the eastern 
United States. Glob ChangBiol. doi:10.1111/gcb.l2152 (2013). 

18. Gloor, M. et al. Does the disturbance hypothesis explain the biomass increase in 
basin-wide Amazon forest plot data? Global Change Biol. 15, 2418-2430 (2009). 

19. Chambers, J. Q. et al. Hurricane Katrina's carbon foortprint on U.S. Gulf Coast 
Forests. Science 318, 1107 (2007). 

20. Kupfer, J. A., Myers, A. T., McLane, S. E. & Melton, G. N. Patterns of Forest 
Damage in a Southern Mississippi Landscape Caused by Hurricane Katrina. 
Ecosystems 11, 45-60, doi:10.1007/sl0021-007-9106-z (2007). 

21. Piatt, W. J., Doren, R. F. & Armentano, T. V. Effects of Hurricane Andrew on 
stands of slash pine (Pinus elliottii var. densa) in the everglades region of south 
Florida (USA). Plant Ecol. 146, 43-60 (2000). 

22. Slater, H. H., Piatt, W. J., Baker, D. B. & Johnson, H. A. Effects of Hurricane 
Andrew on Damage and Mortality of Trees in Subtropical Hardwood Hammocks 
of Long Pine Key, Everglades National Park, Florida, USA. /. Coast. Res. Special 
Issue No. 21, 197-207 (1995). 

23. Blake, E. S., Rappaport, E. N., Jarrell, J. D. & Landsea, C. W. The Deadliest, 
Costliest, and Most Intense United States Tropical Cyclones from 1851 to 2004 
(And Other Frequently Requested Hurricane Facts). NOAA Technical 
Memorandum NWS TPC-4, (NOAA Tropical Prediction Center, National 
Hurricane Center, Miami, FL, 2005), pp 44. [Available online at http://www.nhc. 
noaa.gov/pdf/NWS-TPC-4.pdf] 

24. USDA Forest Service. Forest inventory and analysis national program: U.S. 
Department of Agriculture, Forest Service database, http://www.fia.fs.fed.us/ 
tools-data. (Date of access: 25/11/2012). 

25. Villarini, G., Vecchi, G. A. & Smith, J. A. U. S. Landfalling and North Atlantic 
Hurricanes: Statistical Modeling of Their Frequencies and Ratios. Mon. Weather 
Rev. 140,44-65, doi:10.1175/mwr-d-ll-00063.1 (2012). 

26. Knutson, T. R. et al. Tropical cyclones and climate change. Nat. Geosci. 3, 
157-163, doi:10.1038/ngeo779 10.1038/NGEO779 (2010). 

27. Oswalt, S. N., Oswalt, C. & Turner, J. Hurricane Katrina impacts on Mississippi 
forests. South J. Appl. For 32(3), 139-141 (2008). 

28. Walker, L. R. Timing of Post-Hurricane Tree Mortality in Puerto-Rico. /. Trop. 
Ecol. 11, 315-320 (1995). 

29. Smith, W. B., Miles, P. D., Vissage, J. S. & Pugh, S. A. Forest Resources of the 
United States, 2002. (North Central Research Station, Forest Service - USDA, 
Minnesota, 2004). 

30. Birdsey, R. A. & Heath, L. S. in Productivity of America's Forests and Climate 
Change. In: Gen. Tech. Rep. RM-GTR-271 (Ed. Joyce, L. A.), pp 56-70 (U.S. 
Department of Agriculture, Forest Service, Rocky Mountain Forest and 
Experiment Station, , 1995). 

31. Wang, W., Qu, J. J., Hao, X., Liu, Y. & Stanturf, J. A. Post -hurricane forest damage 
assessment using satellite remote sensing. Agric. For. Meteorol. 150, 122-132, 
doi:10.1016/j.agrformet.2009.09.009 (2010). 

32. Jacobs, D. Potential Damage to Southern Forest Resources 2004-2005 Hurricanes. 
Compass. (October 2008), pp 7 (Science Delivery Group, Southern Research 
Station, Forest Service, 2008). 

33. Chen, S. S., Zhao, W., Donelan, M. A. & Tolman, H. L. Directional Wind-Wave 
Coupling in Fully Coupled Atmosphere-Wave-Ocean Models: Results from 
CBLAST-Hurricane. /. Atmos. Sci. 70, 3198-3215, doi:10.1175/jas-d-12-0157.1 
(2013). 

34. Scheitlin, K. N., Eisner, J. B., Lewers, S. W., Malmstadt, J. C. & Jagger, T. H. Risk 
assessment of hurricane winds for Eglin air force base in northwestern Florida, 
USA. Theor. Appl. Climatol. 105, 287-296, doi:10.1007/s00704-010-0386-4 
(2011). 

35. Emanuel, K. Increasing destructiveness of tropical cyclones over the past 30 years. 
Nature 436, 686-688, doi:10.1038/nature03906 (2005). 

36. FAO. Global Forest Resources Assessment 2005: Process towards sustainable 
forest management. (Food and Agriculture Organization of the United Nations, 
Rome, Italy, 2006). 

37. Lawrence, D. & Foster, D. changes in forest biomass, litter dynamics and soils 
following shifting cultivation in southern mexico: an overview. Interciencia 27, 
400-408 (2002). 

38. Meeker, J. R., Haley, T. J., Petty, S. D. & Windham, J. W. Forest Health Evaluation 
of Hurricane Katrina Damage on the De Soto National Forest. Report No. 2006- 
02-02, (Southern Region, State and Private Forestry, Alexandria Field Office, , 
Atlanta, GA, 2005). 

39. Houghton, R. A., Hackler, J. L. & Lawrence, K. T. The US carbon budget: 
contributions from land-use change. Sri. Cult. 285, 574-578, doi:10.1126/ 
science.285.5427.574 (1999). 



SCIENTIFIC REPORTS | 4:5197 | DOI: 1 0. 1 038/srep05 1 97 



9 



40. Birdsey, R. A., Plantinga, A. J. & Heath, L. S. Past and prospective carbon storage in 
United States forests. For. Ecol. Manage. 58, 33-40 (1993). 

41. HURDAT. National Weather Service Hurricane Center, http://ftp.nhc.noaa.gov/ 
atcf/archive/ (Date of access: 20/09/2012). 

42. Neumann, C. J. The national hurricane center risk analysis program (HURISK). 
(National Hurricane Center, Coral Gables, Florida, 1991). 

43. LPDAAC. Nadir BRDF-Adjusted Reflectance 16-Day L3 Global 500 m. https:// 
lpdaac.usgs.gov/products/modis_products_table/mcd43a4 (Date of access: 19/09/ 
2012). 

44. Schaaf, C. et al. Collection 005 Change Summary for MODIS BRDF7 Albedo 
(MCD43) Algorithms. http://landweb.nascom.nasa.gov/QA_WWW/forPage/ 
C005_Change_BRDF.pdf (2012). (Date of access: 19/09/2012). 

45. Bowyer, P. & Danson, F. M. Sensitivity of spectral reflectance to variation in live 
fuel moisture content at leaf and canopy level. Remote Sens. Environ. 92, 297-308, 
doi:10.1016/j.rse.2004.05.020 (2004). 

46. Texas Forest Service. Hurricane Rita Timber Damage Assessment Texas Forest 
Service. http://texasforestservice.tamu.edu/uploadedFiles/Sustainable/Hurricane_ 
Rita/ritaassessment.pdf (2005) (Date of access: 24/11/2013). 

47. Texas Forest Service. Hurricane Ike Timber Damage Assessment: Texas Forest 
Service. http://texasforestservice.tamu.edu/uploadedFiles/Sustainable/ 
Hurricanelke/ HurricaneIke_Timber-Damage-Assessment.pdf (2008) (Date of 
access: 24/11/2013). 

48. Hall, P. Using the Bootstrap to Estimate Mean Squared Error and Select 
Smoothing Parameter in Nonparametric Problems. /. Multivar. Anal. 32, 
177-203 (1990). 

49. Blackard, J. et al. Mapping U.S. forest biomass using nationwide forest inventory 
data and moderate resolution information. Remote Sens. Environ. 112, 
1658-1677, doi:10.1016/j.rse.2007.08.021 (2008). 

50. Schlesinger, W. H. Biogeochemistry: An Analysis of Global Change, (Academic 
Press, New York, 1991). 

51. Oswalt, S. N. & Oswalt, C. M. Relationships between common forest metrics and 
realized impacts of Hurricane Katrina on forest resources in Mississippi. For. Ecol. 
Manage. 255, 1692-1700, doi:10.1016/j.foreco.2007.11.029 (2008). 



Acknowledgments 



This research, part of the LandCarbon Project, was funded by the Land Change Science 
Program of the U.S. Geological Survey. The authors greatly thank Shengli Huang, Bhaskar 
Ramachandran, and Thomas Adamson for internally reviewing the manuscript. We greatly 
appreciate Prof. Emery Boose from Harvard University, who was consulted on the 
meteorological model of HURRECON. Any use of trade, firm, or product names is for 
descriptive purposes only and does not imply endorsement by the U.S. government. Work 
by SIGT Inc. performed under USGS contract G10PC00044. 

Author contributions 

S.L. and D.D. designed the research. D.D. performed research and wrote the manuscript 
under the guidance of S.L. J.O. contributed on collecting data and writing up the 
manuscript. All authors reviewed the manuscript. 

Additional information 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Dahal, D., Liu, S. & Oeding, J. The Carbon Cycle and Hurricanes in 
the United States between 1900 and 2011. Set Rep. 4, 5197; DOI:10.1038/srep05197 (2014). 

I This work is licensed under a Creative Commons Attribution 3.0 Unported License. 
I The images in this article are included in the article's Creative Commons license, 
unless indicated otherwise in the image credit; if the image is not included under 
the Creative Commons license, users will need to obtain permission from the license 
holder in order to reproduce the image. To view a copy of this license, visit 
http://creativecommons.Org/licenses/by/3.0/ 



SCIENTIFIC REPORTS | 4:5197 | DOI: 1 0. 1 038/srep05 1 97 



10 



